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Comment on "Universal 
Correlated systems" 

In the "universal" probability density function 
(PDF) of various global measures x from correlated equi- 
librium and non-equilibrium physics-inspired models was 
shown (figure 2 of |^ ) to be in reasonable agreement with 
n(y) = K(e}>~'^^ Y , where y = h(x — s). This Bramwell- 
Holdsworth Pinton (BHP) PDF with shape parameter 
a = 7r/2, was contrasted with the Fisher-Tippett-Gumbel 
(FTG) PDF where a = 1, one of 3 possible asymptotes 
for the maxima of uncorrelated random variables (r.v.'s) 
H. In the BHP PDF was also seen to describe the 
distribution of maxima ^^^j. from L trials of length N 
vectors ^ of correlated r.v.'s, generated by ^ = Mx. The 
elements of x were i.i.d. exponentially distributed, and 
M was "an N x N matrix with random but fixed ele- 
ments" . The applicability of the BHP PDF to all the con- 
sidered processes in was inferred to be due to "finite 
size, strong correlations (our italics) and self-similarity" . 
In this Comment we show that this inference is not sup- 
ported by the extremal behaviour of £,max' because 1) 
while a suitable choice of N ensures that ^^ax f^'Hs on 
the BHP curve, other equally arbitrary values of N give 
maxima which don't; and 2) ^ is in fact not strongly 
correlated, and the PDF of ^^ax change little for decor- 
related ^. We thus suggest that a differs from its FTG 
asymptote of 1 here largely due to slow convergence of 
maxima of the nearly Gaussian ^. 

In "random" implies "i.i.d Gaussian, mean 0, 


Fluctuations in effective degrees of freedom due to correlation 
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but fixed" implies M has random elements held con- 
stant between trials. ^ is correlated but the fact that 
M and x statistically independent means that for 
the simplest case of = 2, (= MuXi + M12X2) 
has PDF P(ei) = / dAhidMi2dxidx2S{ii - Mnxi - 
Mi2X2)P{Mii)P{Mi2)P{xi)P{X2)- Taking exponential 
P{x) = Ae^'^^ as in fit] we generalise to N x N matrix M 


so p(o = ^ n m 
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We then have 


is the generating 

lik) = 

f|e^'/2fc^(l - erf(A/^/2k)). 

The maxima of ^ are dominated by the tail behaviour 
of P{^), which from a small k expansion of / is seen 
to be Gaussian ^ (X/^/N) exp —^"^ X"^ /4N . The maxima 
of 5 will thus be in the basin of attraction of the FTG 
distribution. The asymptotic shape parameter a corre- 
sponding to L, TV — > 00 cannot be changed by short range 
dependence , but one expects corrections to a for finite 
N, L because of the logarithmically slow convergence of 
the maxima of Gaussians, and, sometimes, a change in 


Figure 1 shows the PDF of ^4^^ for L = 10^ trials 
in 2 cases: TV = 40 and 250, and for comparison the 
PDF of Xmax = max(x^), when N = 40, to illustrate the 
expected rapid convergence of Xmax the FTG curve. 
Convergence of £,f^ax however slows with increasing N. 
By N = 40, convergence of £,^ax noticeably slower 
than Xmax and the PDF follows the BHP curve; while 
for TV = 250 even slower convergence in Xmax gives ob- 
vious deviation from the BHP form resembling that seen 
in figure 2 of [Q . For TV = 40 the autocorrelation func- 
tion of X (not shown) drops fast by the first lag. Weak 
(less than ±0.2) long ranged oscillations are seen at all 
lags. Despite this, decorrelating by shuffling the L x N 
matrix before taking the maximum of ^ changes the PDF 
of its maxima in a relatively minor way (see Figure 1), 
confirming that correlation cannot explain most of the 
change from the FTG PDF here. Our result implies that, 
even though subsequent results may show that the BHP 
curve can result from strong correlation, it need not. 

NWW thanks Mervyn Freeman, Clare Watt, Pat 
Espy and David Stephenson for valuable suggestions. 
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FIG. \. L = W. Middle plot shows PDF of ^^^^ (cir- 
cles) compared with PDF of maxima of decorrelated process 
(crosses), TV — 40. Top plot (stars) is PDF of maxima Xmax 
of exponentially distributed r.v's, TV = 40. Lowest plot (cir- 
cles) shows ^maa; in TV = 250 case. Also shown on each 
plot are the BHP curve (solid line), and the a = 1 FTG 
plot (dashed line). Displacements of 10~^'^ and 10~^ for the 
middle and lower plots, respectively, were applied for clarity. 
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